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Abstract. An extreme-mass-ratio burst (EMRB) is a gravitational wave signal emit- 
ted when a compact object passes through periapsis on a highly eccentric orbit about a 
much more massive body, in our case a stellar mass object about the 4.3 1 x 1O 6 M mas- 
sive black hole (MBH) in the Galactic Centre. We investigate how EMRBs could con- 
strain the parameters of the Galaxy's MBH. EMRBs should be detectable if the periap- 
sis is r p < 65r g for a pt — 1OM orbiting object, where r g = GM./c 2 is the gravitational 
radius. The signal-to-noise ratio p scales like log(p) -2.7 log(r p /rg)+log( / L//M )+4.9. 
For periapses smaller than ~ 10r g , EMRBs can be informative, providing good con- 
straints on both the MBH's mass and spin. 



1. Introduction 



We believe that most galactic nuclei have harboured a m assive black hole (MBH) during 



their evolution (Lynden-Bell & Rees 1971; Rees 1984). Observations show there are 



corre l ations between the MBHs' masses and their h ost galaxies' properties ( Kormendy & Richstone 
11995b iMagorrian et al.lll998l: iGraham et al.ll201l|) . These suggest coeval evolution of 



MBH and galaxy (Pengl 120071 : Ijahnke & Maccioll201 ll) . Since they share a common 



history; one can inform us of the other. 

The best opportunity to study MBHs comes from the compact object in our own 
galactic centre (GC), coincident with Sagittarius A* (Sgr A*). This is an MBH of mass 
M. = 4.31 x 10 6 M G at a distance of only R = 8.33 kpc dGillessen et al.ll2009h . 

According to the no-hair theorem , the MBH should be described completely by its 
mass M. and (dimensionless) spin a„ (IChandrasekhad [1998 ) . The spin is related to the 
BH's angular momentum / by a* = cJ/GMj. As we have a good estimate of the mass, 
to gain a complete description we have only to measure the spin. 



An MBH accumulates mass and angular momentum through accretion (IVolonteri 

120101) . There are several possible accretion mechanisms, each leaving its own imprint 
on the final spin. Measuring the spin of the MBH shall help us understand the relative 
importance of these processes, and gain insight into the Galaxy's past. 

An exciting means of inferring information about the MBH is through gravi- 
tati onal waves (GWs) emitted wh en stellar mass compact objects (COs) pass close 
by dSafhyaprakash & Schutzl l2009). A space -borne detector can detect GWs in the fre- 
quency range of interest for these encounters (iDanzmann & Rudige i]|2003l ; lAmaro-Seoane et al. 
|2012|) . The w aveforms generated when COs inspiral towards an MBH have been 
much studied (Glampedakis 2005). These systems are typically formed following 



two-body encounters so the initial orbits are highly eccentric; a burst of radiation is 
emitted during each periapse passage. These are extreme mass-ratio bursts (EMRBs; 
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Rubbo, Holley-Bockelmann, & These orbits can evolve, becoming more 



circular, and then emitting continuously i n the detector frequency ra nge. These signals 
are extreme mass-ratio inspirals (EMRIs; lAmaro-Seoane et al.ll2007l) . 

We investigate the bursts from the initial highly-eccentric orbits, under the simpli- 
fying assumption that these orbits are marginally bound, or parabolic, since highly ec- 



centric orbits appear almost indistinguishable from an approp riate parabolic orbitfl The 



event rate for detectable EMRBs has been estimated as 1 yr 1 dHopman. Freitag. & Larson 



2007). Even if only a single burst is detected, this should still give an unparalleled probe 
of the spacetime of the GC. What can be inferred depends upon the orbit, which we in- 
vestigate here. 

We use the classic Laser Interferometer Space Antenna (LISA) design for this 
work. It is hoped that any future detector shall have comparable sensitivity to LISA, 
and that studies using this design shall be a sensible benchmark. 

2. Numerical kludge waveforms 

For given angular momenta, and initial starting position, we can calculate the geodesic 
trajectory in a Kerr background. The orbiting body is assumed to follow this track 
exactly; we ignore evolution due to the radiation, which is negligible for EMRBs. 
From this trajectory we cal culate the waveform using a semirelativistic approxima- 



tion (RuffinijfcSasaki||l981): we assume the particle moves along a geodesic, but ra- 



diates as in flat spacetime. This is known as a numerical kludge (NK), an d well ap 



proximates results computed by more accurate methods (Bab ak et al.1 12007). The use 
of the geodesic ensures the correct frequency components appear in the waveform, but 
the flat-space wave generation means they do not have precisely the correct amplitudes. 

NK approximations aim to encapsulate the main characteristics of a waveform by 
using the exact particle trajectory (ignoring inaccuracies from radiative effects and from 
the particle's self-force), whilst saving on computational time by using approximate 
waveform generation. 

We build an equivalent flat-space trajector y from the Kerr geodesic. This is done 
by identifying the Boyer-Lindquist coordinates ( Boyer & Lindquislll 19671) with a set of 



flat-space coordinates. We use s pherical polars so {^bl, #bl> 0bl) - > {r S ph, #sph, 0sph) 
dGair. Kennefick. & Larsonll2005h . Oblate-spheroidal coordinates yield similar results. 

With a flat-space trajector y, we may use a flat-space wave generation formu la: the 
quadrupole-octupole formula (Bekensteinl ll973l ; iPressI [l977l : lYunes et alj |2008). This 



is correct for a slowly moving source, and is the familiar quadrupole formula, derived 
from linearized theory, plus the next order terms. 



3. Waveforms and detectability 
3.1. Model parameters 

The waveform depends on the properties of the MBH; the CO and its orbit, and the 
detector. We assume the position of the detector is known, and the MB H is coincident 



with the radio source of Sgr A*, which is within 20r g of the MBH (IReid et ajj 2003; 



'Here "parabolic" and "eccentricity" refer to the energy of the geodesic and not to its geometric shape. 
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Doeleman etaDl2008h. We use the J2000.0 coordinates, which are determined to high 
accuracy dReid et al.l ll999: Yusef-Zadeh_etal|U999|). The parameters left to infer are: 
(1) The MBH 's mass M. ; This is well constrained b y the observation of stellar orbits 
about Sgr A* dGhez et al.ll2008l : lGillessen et al.ll2009h . the best estimate is M, = (4.3 1 ± 
0.36) x 10 6 M o . (2) The spin parameter a t . (3, 4) The orientation angles for the black 
hole spin ©k and Ok. (5) The ratio of the GC distance and the compact object mass 
£ = Ro/n- This scales the amplitude of the waveform. Bursts do not undergo orbital 
evolution, hence we cannot break the degeneracy between Rp and u. The distan ce is 
constrained by stellar orbits to be R = 8.33 ± 0.35 kpc dGillessen et al.l l2009h . (6, 
7) The angular momentum of the CO, parametrized by the magnitude at infinity 
and the orbital inclination i. (8-10) Coordinates specifying the trajectory. We use the 
angular phases at periapse, tp v and 9 p , as well as the time of periapse ? p . 



3.2. Signal-to-noise ratio 

The detectability of a burst depends upon its SNR p. The amplitude of the wave- 
form is proportional to the CO mass; we work in terms of a mass-normalised SNR 
p = (ji/M Q )~ l p. We considered a range of orbits. The spin of the MBH and the or- 
bital inclination were randomly chosen, and the periapse distance was set so that the 
distribution would be uniform in log-space (down to the inner-most orbit). For each set 
of the extrinsic parameters, the periapse positions, orientation of the MBH, and orbital 
position of the detector were varied: we used five combinations of these intrinsic pa- 
rameters (each drawn from a separate uniform distribution), taking the mean of lnp for 
each set. 

The correlation between the periapse radius and SNR is shown in fig. CD The shape 
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Figure 1. Mass-normalised SNR as a function of periapse radius. The points are 
the values averaging over each set of intrinsic parameters. The best fit line is log(p) = 
-2.69 log(r p /r g ) + 4.88 for orbits with r p > 13.0r g . 



is predominantly determined by the noise curve. The change in the trend reflects the 
transition from approximately power law behaviour to the bucket of the noise curve. We 
fit a fiducial power law to orbits with a characteristic frequency of /* - ^jGM,/r p < 
1 X 10 -3 Hz, to avoid spilling into the bucket. Changing the cut-off within a plausible 
region alters the fit coefficients by ~ 0. 1 . 

Setting a threshold of p = 10, a 1M Q (10M o ) CO is detectable if r p < 27 r g (65r g ). 
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4. Parameter estimation & results 



Having detected a GW signal, we are interested in what we can learn about the source. 
We performed Markov chain M onte Carlo (M CMC) simulations to characterise the pos- 
terior probability distribution (Mac Kayll2003l . chapter 29). Waveforms were computed 
for a range of orbits. In each case the MBH has the standard mass and position. The CO 
was chosen to be 10M o , as the most promising candidates for EMRBs would be stel- 
lar mass black holes: they are massive and hence produce higher SNR bursts, they are 
more likely to be on close orbits as a consequence of mass segregation, and they cannot 
be tidally disrupted. Orbits were chosen with periapses uniformly distributed in log- 
arithmic space (to the inner-most orbit). The other parameters were chosen randomly 
from appropriate uniform distributions. 

It is possible to place good constraints from the closest orbits, although there is 
significant correlation between parameters. The standard deviation o~ of the recovered 
posteriors are shown in fig.|2]for the mass and spin. Filled circles are used for converged 
MCMC runs. Open circles are for those yet to converge; widths should be accurate to 
a factor of a few. These results do not incorporate any priors (save to keep them within 
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Figure 2. Distribution widths as functions of periapse r p . Filled circles are used 
for converged runs; open circles for unconverged runs. The dotted line is the current 
uncertainty for M.; the dashed line the standard deviation for an uninformative prior, 
and the solid line the total prior range. 

realistic ranges). Therefore, the resulting distributions characterise what we could learn 
from EMRB's alone. 



4.1. Scientific potential 

The current uncertainty in the mass is o~m. - 0.36 x 10 6 M o (~ 8%). It appears that 
orbits of a /u = 10M o CO with periapses r p < 13r g should be able to match this. 
Accuracy of 1% is possible if r p < 8r g . 

The spin is less well constrained. To obtain an uncertainty for the magnitude of 
0.1, comparable to that achieved in X-ray measurements of active galactic nuclei, the 
periapsis needs to be r p < llr g . As the spin encodes information of the merger and 
accretion history, this could illuminate the MBH's formation. 

We have no a priori knowledge about the CO or its orbit, so anything we learn 
would be new. However, this is not particularly useful information, unless we observe 
multiple bursts, and can start to build up statistics for the dynamics of the GC. Using 
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current observations for the distance to the GC, which could be further improved by the 
mass measurement from the EMRB, it is possible to infer a value for the mass /j, from 
f . This could inform us of the nature of the object (BH, NS or WD) and be a useful 
consistency check. A small value of indicating a massive CO, would be unambiguous 
evidence for the existence of a stellar mass black hole. 
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